Evidence of anomalous dispersion of the generalized sound velocity in glasses 
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The dynamic structure factor, S{Q,ij), of vitreous silica, has been measured by inelastic X-ray 
scattering in the exchanged wavevector (Q) region Q=4-^16.5 nm~* and up to energies huj— 115 
meV in the Stokes side. The unprecedented statistical accuracy in such an extended energy range 
allows to accurately determine the longitudinal current spectra, and the energies of the vibrational 
excitations. The simultaneous observation of two excitations in the acoustic region, and the persis- 
tence of propagating sound waves up to Q values comparable with the (pseudo-)Brillouin zone edge, 
allow to observe a positive dispersion in the generalized sound velocity that, around (5~5 nm~^, 
varies from ~6500 to ~9000 m/s: this phenomenon was never experimentally observed in a glass. 

PACS numbers: 63.50.+X, 61.10.Eq, 61.43.Fs 



The stimulating evidence of the existence of short 
wavelength phonon-like excitations sustained by non pe- 
riodic solids, as demonstrated by recent Inelastic X-ray 
Scattering (IXS) studies I*], has motivated a massive 
number of investigations. This notwithstanding, a gen- 
eral consensus on the characteristics of observed excita- 
tions is still lacking 0]. Among several specific ques- 
tions waiting for an answer, one of the most important 
regards the existence and the features of the dispersion 
relations in the whole first pseudo-Brillouin Zone, and 
their connection with the well known universal anoma- 
lies of glasses (low temperature transport and thermody- 
namic properties Q, excess in the density of states 0, 
etc.). 

In the specific case of vitreous silica, the archetype 
of the strong glasses, up to now the experimental stud- 
ies were limited to i) the low Q region ((5<4 nm^^) - 
investigated by IXS- where the spectra are described by 
a single excitation with energy position {hfl{Q)) linearly 
dispersing with Q 0, Ei 0, and ii) the high Q region 
{Q>15 nm^^) where Inelastic Neutron Scattering (INS) 
data (confined to the high Q and low uj portion of the 
S{Q, w) due to the kinematic restrictions) show again 
spectra dominated by a single excitation, but having now 
a Q independent energy This scenario led to the 

conclusion that the longitudinal acoustic branch in vitre- 
ous silica -dispersing with a sound velocity ww6400 m/s 
at low Q- fiattens at increasing Q giving origin to the Bo- 
son Peak 9]. However, as the cross-over between the low- 
{il{Q)=vQ) and high-Q {^{Q)=^bp) regimes lies in the 
not-explored Q region, it was not possible to establish a 
firm conclusion on the shape of the dispersion relation 
of this mode. A hint on this issue came from Molecular 
Dynamic (MD) simulations. Surprisingly, the MD works 
showed the simultaneous presence of two excitations in 
the longitudinal current spectra at large enough Q {Q>5 



nm-i) [i3,IilElIP|- MD indicated that the excitation 
detected at low Q by IXS corresponds to the Longitudinal 
Acoustic (LA) mode, while that detected at high Q by 
INS is associated to the spilling of the Transverse Acous- 
tic (TA) mode in the longitudinal spectra and not to the 
longitudinal branch 13]. The presence of the signature 
of the transverse dynamics in the longitudinal spectra is 
not surprising. Indeed, in a disordered medium, the po- 
larization character of the modes becomes ill defined at 
short wavelengths -when the sound waves no longer see 
the system as a continuum- and this gives rise to modes 
with mixed polarization 

The MD works [H E [13 also indicated the pres- 
ence, in vitreous silica, of another important feature of 
the LA branch: a positive dispersion of the general- 
ized sound velocity, v{Q)=Q{Q)/Q. A positive disper- 
sion of v(Q) is usually observed in viscoelastic liquids 
where a step-like change of v{Q) is found when the con- 
dition n{Q)Ta{Q,T)^l is fulfilled, being Ta{Q,T) the 
structural (a) relaxation time. This condition marks the 
transition from the low frequency (and low v) viscous 
behavior to the high frequency (and high v) elastic be- 
havior of the liquid. In glasses, being Tq(Q,T) basically 
infinite, the MD observed positive dispersion cannot be 
associated to the structural relaxation process but to a 
further relaxation process named instantaneous or micro- 
scopic. This process, already observed in the MD simula- 
tion of monatomic glasses L5l and, experimentally, in liq- 
uid metals above melting [ig, has been so far explained 
in term of structural disorder jl^j . No experimental evi- 
dences for such a positive dispersion in glasses have been 
reported yet. 

In this Letter we present an IXS study of the dynamic 
structure factor of vitreous silica in the intermediate Q 
region (4<Q<16.5 nm^^). Using also existing IXS data 
at (5<4 nm""'^, we draw a self-consistent picture of the 
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whole pattern of excitations in the first pseudo-Brillouin 
Zone (p-BZ) of this prototype glass. Specifically, we find 
that i) the LA branch spans energies as high as 80 meV, 
confirmingand extending the validity of previous IXS 
studies lEim limited to (5<4 nm"^ and to E<20 meV; 
ii) above Q=4 nm~^ the spectra show the simultaneous 
presence of two excitations, the high energy one associ- 
ated to the LA modes, and the low energy one assigned 
to the TA excitation. The latter excitation appears in 
the longitudinal spectra due to the mixing phenomenon 
|l4|. and its experimental observation confirms previous 
MD results lldflllll^ and parallels similar findings 
in glassy glycerol . iii) Finally, and most importantly 
the clear identification of longitudinal acoustic branch 
allows to observe experimentally the existence in glasses 
of a positive dispersion of the generalized longitudinal 
sound velocity. The transition is found at Q~5 nm^^, 
where v undergoes a change from «6400 to «9000 m/s). 

The measurements were performed in two different ses- 
sions at the inelastic X ray scattering beamline ID 16 of 
the European Synchrotron Radiation Facility (ESRF). 
The experimental setup consists of a backscattering op- 
erating monochromator and of a 6.5 m analyzer arm host- 
ing a five analyzers bench with a constant angular offset 
between each analyzer. The incident energy was 17794 
eV, corresponding to the Si(9,9,9) reflection. In this con- 
figuration the overall (FWHM) energy and Q resolutions 
were SE = 3.0 meV and 5Q = 0.27 nm~^, respectively. 
The exchanged wavevector is selected by rotating the an- 
alyzer arm. Each energy scan was performed at constant 
momentum transfer by varying the relative temperatures 
of the monochromator and the analyzer crystals. Further 
details on the beamline are reported elsewhere 0| . The 
sample of vitreous sihca (SUPRASIL), heated at T=1270 
K, was studied in a Q range not covered by previous mea- 
surements, i. e. at (5=4, 5, 6 in a first run and from 6.5 to 
16.5 nm~^ in step of 2.5 nm"^ in a second run. In order 
to follow the evolution of the inelastic signal, in the sec- 
ond run we performed wide energy scans up to 115 meV. 
To achieve the necessary statistical accuracy the spectra 
were collected on the Stokes side only (and on a small 
portion of the anti-Stokes to cover most of the central 
peak): the total integration time was 30 minutes/point, 
i. e. Ril20 hours for the collection of a whole spectrum. 

As an example, the spectra obtained for three selected 
Q values are shown in Fig. 1. The IXS signal (open 
circles), proportional to the dynamic structure factor, 
S{Q,iu), convoluted with the experimentally determined 
instrumental resolution function, R{uj), is reported in log 
scale together with the elastic (dotted line, coincident 
with R{u})) and inelastic (dashed line) contributions as 
derived from the fit (see below). Directly from the IXS 
data -at least for the lowest two reported Q values- one 
can infer the presence of two shoulders, which are not 
present in the smoothly decreasing resolution function. 
These shoulders, indicated in the figure by arrows, point 
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FIG. 1: IXS spectra of vitreous silica at r=1270 K for se- 
lected values of the exchanged wavevector in the region of 
interest (open circles with error bars). The total fitting result 
with two excitations (full line) is also reported, along with the 
elastic (dotted line) and inelastic (dashed lines) contributions. 
The arrows indicate the approximate energy location of the 
excitations "bumps". 



out the presence of two excitations. To get quantita- 
tive information on the energies of the excitations and 
on their Q-depcndences, the data have been fitted by the 
convolution of R{uj) with a model function for S{Q,lu) 
weighted by the detailed balance function. The model 
function is the sum of an elastic and an inelastic con- 
tribution, S{Q,U!) = Sei{Q,oj) + Sinei{Q,^)\ the former, 
accounting for the frozen structural relaxation, is repre- 
sented by a delta function, the latter has been tentatively 
described by one or two excitation(s) model through, re- 
spectively, one or two Damped Harmonic Oscillator func- 
tion(s). This model function has been chosen only to ex- 
tract information on the position of the peaks. Indeed, a 
full description of the spectral shape in the investigated Q 
region would need a detailed generalized hydrodynamic 
model, where the transverse and the longitudinal vari- 
ables couple each other due to the broken ergodicity in 
the glass. The use of this detailed model to describe the 
spectra goes beyond the aim of the present Letter. 

All the measured IXS spectra have been fitted using 
both one and two inelastic features, and the F-test has 
been applied to the results in order to establish whether 
the second excitation is statistically necessary to describe 
the spectra: this is always the case for Q values larger 
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FIG. 2: The current spectrum of the inelastic part of the IXS 
signal (the elastic contribution to the fit is subtracted by the 
experimental data) at (5=6.5 nm~^ (open circles) is reported 
together with its ±(j error bars and the one- (dot-dashed) and 
two- (dashed) excitations fits. The inset shows the current 
spectra of the whole IXS signal, without subtraction of the 
elastic contribution, and the total fit to the spectrum (dashed 
line) including the elastic and inelastic (two- excitations) part. 

than 4 nm"-'^. The statistical significance of the two exci- 
tations fit of the spectrum at (3=4 nm~^ is at the thresh- 
old. Previous IXS spectra at lower (5's have been re- 
analyzed and no evidence for the second peak has been 
found. 

It is worth to note that, due to the weakness of the 
inelastic signal with respect to the intense elastic peak in 
this sample, the difference between the one- and two- 
excitations fits could not be appreciated from Fig. 1. 
However, beyond the quantitative statistical data anal- 
ysis, the presence of two excitations clearly arise look- 
ing at the longitudinal current spectra, Cl{Q,^) — 
{uj'^/Q'^)S{Q,llj), as it is evident in Fig. 2. As an ex- 
ample, in this figure we report for (3=6.5 nm~^ the dif- 
ference, multiplied by w^, between the IXS data and the 
Sei{Q,i^) determined by the fit (open circles with error 
bars). The broken lines represent uj'^Si„ei{Q,uj) obtained 
by the previously described fit of the measured spectra 
with one (dash-dotted) or two (dashed) excitations. As 
can be easily seen, one single excitation does not account 
even qualitatively for the measured spectrum. Even if 
the agreement between the data and the two excitations 
fit is not perfect, pointing out the need for a more com- 
plete generalized hydrodynamic model, it is clear that a 
two-excitations description of the spectra is compulsory 
in this Q region. 

To show that the need for a two excitation represen- 
tation of the data -as demonstrate in the main panel of 
Fig. 2- is not an artifact due to the subtraction of the 
elastic line, the " experimental" current spectra (raw IXS 
data multiplied by w^) is reported in the inset of Fig. 2. 
The presence of the (almost Lorentzian) resolution func- 



tion manifests itself as a third bump at low energies (wl.S 
meV), and as a long (almost w-independent) tail. Beside 
these two resolution features, the presence of two excita- 
tions is evident also in the inset of Fig. 2, i. e. in the raw 
data. 
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FIG. 3: A) Dispersion relation of vitreous silica in the first 
p-BZ. Full symbols are experimental data (maxima of the 
two inelastic contribution to the current spectra) representing 
the excitation energies {hi}{Q)) of the longitudinal (diamond) 
and transverse (circle) branches. Data from (3=4 nm~^ are 
from present work, below Q—4 nm~^ are from Ref.s 0, 
Q (triangles). The open symbols are the main maxima of 
the simulated L (open diamond) and T (open circle) current 
spectra from Ref. fl3l . Dash-dotted and dashed lines are the 
extrapolation of the low frequency L and T sound velocity, 
respectively. B) Generalized sound velocity derived from the 
data of panel A as v{Q)=Q,{Q)/Q. Symbols as in A, the MD 
data are not reported. 

The excitation energies obtained by the fit and the cor- 
responding apparent sound velocities {v{Q)=Q{Q)/Q) 
are reported in Fig.s 3A and 3B respectively. Data from 
previous IXS measurements in vitreous-silica 0, S 01' 
performed at low Q are also reported. The corresponding 
dispersion curves obtained from numerical studies [l^ as 
the maxima of the longitudinal and transverse currents 
are plotted in the same figure (open symbols). It is ev- 
ident from Fig. 3 that the two excitations found by the 
fit analysis of the experimental data match with the LA 
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and TA branches derived from the MD simulations. In 
particular, the lower energy excitation in the IXS data 
is in very good agreement with the TA branch, and is 
then associated with the spilling of the TA excitation in 
the measured (longitudinal) spectra due to polarization 
mixing 0| that, in this glass, becomes effective for Q's 
larger than RiS nm~^. 

The excitation at higher energy obtained from the fit 
of the measured spectra is associated to the LA mode: 
the ri((5) values obtained here are the extension at high 
Q's of the LA dispersion relation previously determined 
by IXS and are also in agreement with the trend 

of the LA branch derived from MD simulations. This 
branch shows a clear bend-up at w5 nm~^, i. e. the gen- 
eralized sound velocity undergoes a speed up from its 
low Q value (t;«6400 m/s) to «9000 m/s [l^- The ex- 
istence of a positive dispersion of the generalized sound 
velocity could be interpreted as the signature of an un- 
derlying relaxation process. This possibility, as well as 
the origin of this process -that cannot be identified with 
the structural (a) relaxation process, basically frozen in 
the glassy phase- is an open problem. In Ref. it has 
been suggested that the origin of this relaxation process 
could be retrieved in the nature of the spatial pattern 
of the vibrational eigenmodes of a disordered structure, 
rather than in dynamical effects as, for example, the an- 
harmonicity of the vibrational dynamics. Further study, 
specifically on the temperature dependence of this phe- 
nomena, are required to asses this point. It is important 
to mention that, by a INS study of the dynamics of v- 
silica at large momentum transfer {Q >20 nm~^), Aral 
et al. [i^l deduced a generalized sound velocity of «9400 
m/s, a value compatible with the one found here in the 
high frequency side of the transition. Finally, at Q larger 
than Ri7-8 nm~^, the LA dispersion relation shows the 
features of the usual acoustic phonons behavior: it starts 
to bend down and reaches a maximum at half of the p-BZ 
edge (wll nm~^ in vitreous silica). 

In conclusion, in the present work we have studied the 
dynamics in the intermediate Q region (from 20% up to 
the first p-BZ edge) in vitreous silica. Wc have shown 
that, beside the usual propagating longitudinal sound 
modes already observed in glasses at "small" Q's (up 
to 20% of the p-BZ edge in i)-silica), a second excitation 
becomes more and more evident in the spectra as Q is 
increased. This second excitation -in the Q range where 
it is visible- is almost non-Q dispersing and is assigned to 
the Transverse Acoustic (TA) branch. The possibility to 
observe simultaneously the TA and LA excitations has 
allowed the detailed determination of the shape of the 
LA branch, which shows three distinct regions: i) At Q 
below 4 nm^^ it linearly disperses with a sound velocity 
wRi6400 m/s, consistently with light scattering measure- 
ments [i^. ii) with increasing Q, it does not show ten- 
dency to saturation 0, but, on the contrary, it reveals 
a speed-up, and it reaches a generalized sound velocity 



of «9000 m/s. iii) At Q higher than 7 nm"-'^, it starts 
to follow a crystal-like behaviour with a slow down and 
a maximum around half of the p-BZ. The observation in 
the LA branch of vitreous silica of a positive dispersion 
of the generalized sound velocity provides an important 
missing tile in the phenomenology of disordered induced 
features in glasses. As a speculation, a possible explana- 
tion for such a positive dispersion can be searched in the 
presence of a very fast "relaxation" process associated 
to the structural disorder, as suggested in Ref. or 
in the interaction of the LA branch with other modes, a 
phenomenon similar to the symmetry avoided crossing in 
crystals. 
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